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metries for lepton-pair production in association with a large transverse momentum jet at 
^ large hadron colliders. We find that the leading order results are essentially confirmed. Al- 



We consider the next-to-leading order corrections, 0(a s ), to forward-backward charge asym- 



though experimentally challenging and in practice with large backgrounds, these observables 
could provide a new determination of the weak mixing angle sin 2 9^ t (M^) with a statistical 
precision for each lepton flavour of ~ 10~ 3 (7 x 10~ 3 ) at LHC (Tevatron), and if b jets are 
identified, of the b quark Z asymmetry Ap B with a statistical precision of ~ 2x 10~ 3 (4x 10~ 2 ) 
at LHC (Tevatron). 
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INTRODUCTION 



The large cross sections for gauge boson production at the Fermilab Tevatron and the 



CERN Large Hadron Collider (LH 



C) might give a chance to determine the electroweak 



parameters with high precision 

HQ. 

In practice, the experimental challenge is very de- 
manding, but at any rate pursuing such measurements will help to disentangle the strong 
physics contributing to the different processes. In this paper we calculate the forward- 
backward charge asymmetries of lepton pairs in events with a large transverse momentum 
jet pp — > Z, 7* + j — > e~e + + j 2 at next-to-leading order (NLO), 0(a s ) corrections. We 
make use of the Monte Carlo program MCFM v4.1 3], which includes the necessary pro- 



cesses at this order, and of ALPGEN 4]. The particularly interesting case of a final b jet 
is discussed in detail. We find that the leading order (LO) predictions |5| are essentially 
confirmed. 

Electron-positron pair production pp* — > Z, 7* — > e~e + has a large cross section at large 
hadron colliders, and as it is sensitive to the presence of vector and axial-vector fermion 
couplings to neutral gauge bosons, in principle allows for their precise measurement. A prime 
example is the determination of the effective weak mixing angle sin 2 O 1 ^, that enters in their 
definition, the optimum observable being the forward-backward charge asymmetry of the 
lepton pairs Afb jslQl- Indeed, the tree level Drell-Yan 3 parton process qq — > Z, 7* — > e~e + 
gives an asymmetric polar angle electron distribution relative to the initial quark, which 
also depends on the lepton pair invariant mass M e - e + 

At the Fermilab Tevatron Run I the Collider Detector at Fermilab (CDF) reported an 
asymmetry A FB at the Z peak of 0.07±0.02 [9|, in agreement with the Standard Model (SM) 
prediction. A new measurement of neutral gauge boson production in pp collisions at the 
upgraded Run II Fermilab Tevatron operated at y/s = 1.96 TeV has been recently presented 
by CDF, giving, at the Z peak, A FB = 0.07 ± 0.03, the statistical error being large because 
only an integrated luminosity of 72 pb _1 has been analysed less than a tenth of the 
luminosity collected so far. 4 A fit to these data where the quark and electron couplings to 



2 Throughout the paper we will explicitly refer to the e e + (+j) decay channel. The same analysis applies 
to + j production with I = /i. 

3 A list of Drell-Yan cross section measurements at Tevatron run I and II is given in Ref . |8j . 

4 In Run I the asymmetry found for the last bin, M e - e + G [300, 600] GeV, deviated from the SM prediction 
what is not confirmed by Run II 0] . 
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the Z boson are expressed as a function of sin 2 O 1 ^ 1 gives sin 2 fljf = 0.2238±0.0040±0.0030, 
where the errors stand for statistics and systematics, respectively. This is far away from the 
estimate of the expected statistical precision ~ 0.0005 to be reached at Run II with an 
integrated luminosity of 10 fb _1 h|. Even if the experiment is well understood [llj and 
the theoretical calculations not ambiguous, it seems very hard to get rid of systematic 
uncertainties to the required level. In the following with a more exclusive process we will 
need a more demanding experimental performance. But we will stay on the very optimistic 
side, emphasizing what we may learn if we were only limited by statistics. At LHC with an 
integrated luminosity of 100 fb _1 the weak mixing angle precision would be hopefully further 
improved by a factor ~ 3. This would be comparable to the current global fit precision, 
0.00016, but for instance a factor ~ 2 better than the effective weak mixing angle precision 

n 

obtained from the bottom forward-backward asymmetry at LEP and SLD, 0.00029 [12J. 

The associated production of a neutral gauge boson V — Z, 7* and a jet j has also a large 
cross section, especially at LHC, and can also allow for a precise determination of sin 2 
This NLO correction to V production is a genuine new process when we require the detection 
of the extra jet. In particular, gluons can be also initial states, and the large gluon content 
of the proton at high energy tends to increase the Vj production cross sections, although 
they stay almost one order of magnitude smaller than the corresponding V cross sections. In 
Table ITT1 of the section that collects our numerical results we gather the different LO and NLO 
contributions to V(— > e~e + )j production at Tevatron, to be compared with the inclusive 
LO and NLO V — > e~e + cross section for the same cuts, 127 and 158 pb, respectively. At 
LHC we find for e~e + production 685 and 745 pb, respectively, to be compared to the e~e + j 
cross section, 53 and 57 pb, in Table |H] below. All the calculations throughout the paper 
have been performed with MCFM v4.1, and with ALPGEN when necessary. They provide 
a good description of these processes at hadron colliders. For instance, the prediction at 
Tevatron for the ratio of the inclusive cross section for pp — > V b to pp — > Vj production 
is, according to the results in Table El below, 0.020 to NLO 5 (0.0096 to LO [5]). 6 This has 
to be compared to the recent measurement of this ratio with the DO detector 0.023 ± 0.005 
13], obtained with a similar, but not identical, set of cuts. 



5 We neglect in this estimate the small fraction of events where b and b combine into the same jet. 

6 In apparent agreement with the NLO prediction 0.018 ±0.004 by J.M. Campbell and Willenbrock quoted 

inRef. 0. 
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As pointed out in Ref. \^ the forward-backward charge asymmetry of the lepton pairs can 
be measured in neutral gauge boson production with an accompanying jet either relative to 
a direction fixed by the initial state A FB as in the inclusive neutral gauge boson production 
(Drell-Yan case), or relative to the final jet direction A FB . The former is adapted to obtain 
the asymmetry from qq events, and the latter from g { q ones. Both asymmetries give similar 
precision for sin 2 at LHC but not at Tevatron, where the precision for A FB is almost 
one order of magnitude higher. However, in principle A FB also allows for the measurement 
of flavour asymmetries. Thus, if we require the final jet to be a b quark, we can make a 
new measurement of A FB . This is especially interesting given its observed deviation at the 
Z pole from the SM prediction, 3 a [2]. However, to approach a similar precision will be 
a very demanding experimental challenge because we have not only to identify the heavy 
quark but to measure its charge. Being very optimistic the corresponding effective weak 
mixing angle precision to be in principle expected at LHC, ~ 10 -3 , is already lower than the 
one reported by LEP and SLD, 2.9 x 10~ 4 , but similar to the difference between the central 



values resulting from A FB at the Z pole and the global fit to all data 



121. 



In the following we study the LHC and the Tevatron potentials in turn. First, we review 
the LO contributions to pp } — > Vj production and introduce the different asymmetries. 
Afterward we discuss the NLO corrections, paying special attention to the case of a final b 
jet. Finally, we present the numerical results and draw our conclusions. 



LO PROCESSES AND FORWARD- BACKWARD CHARGE ASYMMETRIES 



Let us thus compare the processes 

pp ] —>■ Z, 7* —>■ e~e + (1) 

and 

pp -y Z,-y* +j -»• e~e + + j (2) 

at LO and define the different forward-backward asymmetries we are interested in. Fig. ^ 
shows the LO diagram contributing to the Drell-Yan process in Eq. ^ In the absence of 
gluonic radiation the transverse momentum of the exchanged vector boson is zero. Therefore, 
one must, for instance, expect that the direction of the initial quark state and the final e~ 
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FIG. 1: LO qq contribution to the Drell-Yan process in Eq. ^ 

are correlated. When the initial quark line emits gluons, such correlations, although still 
present, tend to diminish because of the transverse momentum p t acquired by the e~e + 
system. If one of those additional gluons is hard enough and is emitted in the central region 
of the detector, it gives rise to an extra jet resulting into the process in Eq. |21 However, 
allowing for an extra jet obliges to consider new subprocess initiated by a gluon and a(n) 
(anti)quark. We show in Fig. |2]the relevant tree level diagrams in the simple case j = b. 
In diagram (b) the decay products of the Z, 7* system know very little about the direction 



9. 

^" b 



\ 



6 — — 



b b 

(a) (b) 



FIG. 2: LO gb contributions to the process in Eq. E] for j = b. 

of the final state b, because of the initial state gluon that also connects to the b fermionic 
line. Instead, in diagram (a), the pt of the b quark exchanged in the s-channel is zero, 
therefore one expects correlations between the final state leptons and the direction of the 
6-jet 7 0. An optimal observable to quantify such correlations for the process of Eq. EJis a 
forward-backward asymmetry: 

F-B 



with 



1 da r° da 

dcos6* B= -dcos6*. 



dcos# J-i dcos6* 



7 Analogous considerations also hold in the general case with j ^ b. 
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One can consider two possible angles: 



cos 9, 



cs 



COS 9n 



2(p e z E e+ -pfE e 



(p e — p e+ ) ■ p> 



(p e ~ + p e+ ) ■ pi ' 

where the four-momenta are measured in the laboratory frame and p? = (0,^,^,0). The 
Collins-Soper angle [3| 9qs is, on average, the angle between e~ and the initial quark 
direction, while 9j is the angle between e~ and the direction opposite to the jet in the e~e + 
rest frame From the previous discussion it should be clear that the former choice is 
adapted to the qq collisions and the latter to the g q ones. 

Different asymmetries can be defined, according to the scheme given in Table HJ A 



Collider 


Asymmetry 


Definition 


PP 


ACS 
^FB 


cos 6 


= cos 


9cs 


pp 


ACS 
^FB 


cos 6 


= cos 


9 x ip^+pr+pii 


PP 


A j 
^FB 


cos 6 


= cos 


9 x \p% +pI~+p p z\ 
3 Pt + +Pt~+pl 


pp 


^FB 


cos 6 


= cos 




PP 


A b 
^FB 


cos 6 


= cos 


9j x (-sign(Qb)) 


PP 


A b 
^FB 


cos 6 


= cos 


9j x (-sign(Q fe )) 



TABLE I: The definitions of the various asymmetries at the pp and pp colliders. 



comment is in order, with respect to the phases appearing in the Table. In pp colliders the 
quark direction is fixed by the rapidity of the jet plus the lepton pair. This implies defining 
cos 9ns with an extra sign factor ^ Pz _ +Pz , +p | I , as in the second line of Table IJ. On the other 
hand, in pp colliders there are produced as many quarks as antiquarks and Ap B vanishes 
unless some difference is made between them. Hence, cos# is defined with an extra sign 
factor — , p = p e +p e+ + p°\ which corresponds to assume that the largest rapidity parton 

Pz 

is a (anti)quark if it is along the (anti)proton direction. This explains the factor in the third 
line of Table HI Finally, because both in pp and pp colliders are produced as many b as b, 
in order to obtain Ap B one must use cos9j multiplied by a +(— ) sign for b (anti)quarks, 
— sign(Q ft ) with Q b the b charge. In practice this means detecting the charge of the produced 
b jet. 
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Such asymmetries have been studied in detail, at LO, in Ref. |5|. What is important to 
realize is that, in order to get reliable predictions, a priori small additional contributions 
must be carefully taken into account. Let us study, in particular, the effect of including 
radiative 0(a s ), NLO, corrections. 



NLO CONTRIBUTIONS 

In the following we discuss all contributions necessary to compute the process in Eq. 121 
at the NLO. The described structure is implemented in the MCFM code, that we used as 
it is in the case j ^ b. However, as already pointed out, the computation of Ap B requires 
disentangling b from b final states. In MCFM this selection is not possible on an event 
by event basis, because b and b contributions are summed up. Therefore, we modified the 
code to take this into account. In addition, we included part of the remaining real NLO 
contribution with the help of ALPGEN. We find convenient to list the different contributions 
in the following, for the case j = b, because this will allow us to discuss their relative size. 
Analogous considerations apply to the general case. 

The virtual contributions are drawn in Fig. El together with the definition of all the 
graphical symbols and conventions used. In particular, we omit drawing explicitly the decay 
of the Z, 7* system, but we always understand Z, 7* — > e~e + , and the blob stands for the sum 
of all possible contributing Feynman diagrams. For the case at hand, it means exchanging 
a virtual gluon in all possible ways in the two diagrams of Fig. El 

b 

» 

9 

q or q 

Z, 7* — > e~e + 




FIG. 3: NLO virtual gb contributions. 

The real contributions are given in Fig. |3J The shorter lines on the right part of the 
drawings means that the corresponding outgoing partons are not seen because they are too 
soft or too collinear to the ingoing or outgoing b quark, therefore also leading to a final state 
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(a) (b) 





(d) 

FIG. 4: NLO real gb (a), q b (b), gg (c), qq (d) contributions. 

formed by an e~e + pair plus a jet containing a b quark. The NLO fully differential cross 
section d<r NLO is given by 



d(T NLO = £ f X dxx f 1 ^ /NLO (;Ei)/ NLO (x2) d ~ 



NLO 



*,j = 9, 9 ( 4 ) 

where /j NLO (a;) is the parton density, evoluted at the NLO, relative to the i th initial state 
particle, carrying a fraction x of the proton (or antiproton) longitudinal momentum 

ff LO {x)=it\x)+a s ft\x) ) (5) 

and do"]j LO are the differential cross sections corresponding to the subprocesses given in Figs. 
E3 and |U computed at the one loop accuracy in QCD 

da^ LO = da§ ) + aMf- (6) 
The separation between virtual and real contributions in do - ,-. 1 is performed in MCFM with 

n 

the help of a Dipole Formalism |15| . 

It is worth discussing the relative size of all the terms appearing in Eq. |3J Setting 
the scale of the hard scattering to fi = Mz, the distribution function fjf\x) intrinsically 
sums up all contributions of the order a k s L k where L = ln(///m&) is the large collinear 
logarithm associated to the fact that fjf\x) describes an exactly collinear gluon splitting 
g — > bb. Therefore, the leading order term 

d*<°> = C dx! C dx 2 fS i \x 1 )fi°\x 2 ) d^ } + (g~b), (7) 



o Jo 
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coming from the diagrams in Fig. |21 contains all possible a k s L k contributions 

oo 

do® = £c fe * k s L k . (8) 

k=l 

Including the diagrams in Figs. 01 Ufa) and Efb) takes into account all corrections order 
a^ k+1 ^L k , because the corresponding subprocesses multiply fl ,(x) in Eq. 0J The corrections 
given in Fig. 0Jc) correspond to the contributions 0(a s ) or 0(a s L), depending whether the 
final state b is or is not collinear to one of the initial state gluons 8 . They correct the leading 
order picture implicit in Eq. [7| of exactly collinearly produced b: now the initial state b can 
acquire a p t . However, as discussed, all collinear contributions are already included in the 
LO process, in particular the k — 1 term in Eq. [SI therefore we are facing an apparent 
double counting. The key for understanding that this is not the case is noticing that other 
contributions are present in Eq. that correspond to the evolution of the parton densities 

«. (fPMfPM + /f(*i)A (1) (* 2 )) d*£> + (g~b). (9) 

Their effect is, among others, subtracting the c\a s L term from the LO contribution in Eq. 
IH1 so that adding the corrections of Fig. Efc), does not imply double counting. The structure 
described so far is implemented in MCFM, that we had to modify in order to disentangle b 
and b production, which is necessary, in our case, for computing A\ B . 

Finally, we computed the pure 0(a s ) corrections in Fig. Hfd) with the help of ALPGEN. 
In the next section we discuss our numerical findings. 



NUMERICAL RESULTS 



We present our numerical results for e e + j and e e + b at LHC and Tevatron in turn. Our 
simulation of the set up at LHC (Tevatron) is as follows 

V\ > 20 GeV, p{ > 50 (30) GeV, 
\r] e ' j \<2.5, AR eJ >0A. 



For muon pairs the main difference would be the 
the cteq611 (cteq6m) parton distributions at LO (NLO) 



pseudorapjdity coverage [13, [l8| . We use 
19] ]. The effect of smearing the 



Note that such collinear or almost collinear configurations contribute to the exclusive process we are 
studying when the final state b is lost in the forward or backward regions of the detector. 



9 



lepton and jet energies has been studied at LO and found to be negligible p, therefore 
we do not include it here. On the other hand, the dominant background processes are 
expected to be the same as for Drell-Yan production, namely, jets misidentified as e 1 * 1 , and 
pp ] — > W + W~j — * e + e~v e v e j . They are understood experimentally, at least at Tevatron, 
and can be considered under control Q|. The LO and NLO production rates at LHC and 
Tevatron are given in Table ITT1 In Figs. 03© we show the corresponding charge asymmetries, 



Contributing 
process 


LHC 
LO NLO 


Tevatron 
LO NLO 


gq Vj(j) 
QQ -> Vj(j) 
qq Vj(j) 
99 Vj(j) 


44.3 
8.4 




3.4 
>3.7 


3.40 
4.61 




4.77 
> 2.76 


Total 


52.7 


57.1 


8.01 


7.53 


9b -> Vb{g) 
99 -> Vb(b) 
fb^VbCq") 
qq - Vb(b) 


1.81 


> 

C 


> 1.81 
.06 


0.038 


> 

C 


> 0.049 
.025 


Total 


1.81 


1.87 


0.038 


0.074 



TABLE II: Estimates for the e~e + j and e~e + b cross sections at LHC (y/s = 14 TeV) and Tevatron 
(i/i = 1.96 TeV) in pb. The jet transverse momenta are required to be larger than 50 (30) GeV 
at LHC (Tevatron) and all pseudorapidities \rj\ smaller than 2.5. The pt of the leptons is larger 
than 20 GeV. The separations in the pseudorapidity-azimuthal angle plane satisfy AR > 0.4 and 
M e - e + is within the range [75, 105] GeV. V means summing over q and q contributions. 

Apg relative to the initial parton and A FB , Ap B relative to the final jet. 

The effect of the 0(a s ) corrections is moderate for Apg and A FB , but sizable for A FB , 
especially at Tevatron. This is basically due to the genuine new higher order process qq — ► 
Vb(b) in Fig. Efd), that tends to wash out the asymmetry which is mainly associated to the 
(a) contribution in Fig. EJ This change is more pronounced at Tevatron energies, where the 
qq content of the (anti)proton is larger (see Table ITTj) . In all cases the asymmetries at NLO 
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FIG. 5: NLO (solid histogram) and LO (points) asymmetries at LHC. 

diminish, except for Ap B at Tevatron where this asymmetry is smaller. 

Near the Z pole, M e - e + ~ Mz, the asymmetries can be approximated by [f| 

A = b(a-sin 2 ^ pt (MD), (10) 

translating then their measurement into a precise determination of sin ^(Ml). In Table 
IIHI we collect the asymmetry estimates, their statistical precision, the cross sections and the 
precision reach Ssii^O 1 ^ of LHC and Tevatron for M e - e + in the range [75, 105] GeV, by 
assuming an integrated Luminosity L of 100 (10) fb^ 1 at LHC (Tevatron). In the Table 
we assumed a 6-tagging efficiency e of 100 %, no contamination u and, in particular, no 
charge misidentification. The statistical precisions 5 A and S sin 2 are proportional to 
e~5 ; and the asymmetries A to 1 — 2cu. Therefore the contamination multiplies 5sm 2 9 ] ^ t 
by (1 — 2uj)~ 1 . Thus, if we only consider semileptonic b decays, implying e ~ 0.1 and 
uj ~ 0, 5 A and ^sin 2 ^^ increase by a factor ~ 3. In practice we must try to maximize 
the quality factor Q = e(l — 2oj) 2 20J. Although the exact value of to can only be inferred 
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FIG. 6: NLO (solid histogram) and LO (points) asymmetries at Tevatron. 

from dedicated experimental studies, it is instructive to have an idea of its typical size. 
This can be achieved by comparing the charge separation 5b, measured at LEP using a jet 
charge technique, with the value 2Q b — — ~ one would get if the b quark could be directly 



observed |2l|. This comparison gives uo ~ 0.3. The statistical precisions given in Table UTT 
are certainly optimistic for systematic errors are also sizable. At any rate approaching the 
quoted precisions will be an experimental challenge. 

A second source of uncertainty, that is not accounted for in Table 1111} is the dependence 
of the results on the chosen set of parton densities. We investigated it by recomputing 
asymmetries and statistical precisions using different parton distribution sets in the classes 



cteq and mrst 22]. By doing so, variations of the asymmetries of the order of 10% can be 
easily observed in the range 75 GeV < M e - e + < 105 GeV, both at Tevatron and LHC, while 
the statistical precisions are not significantly affected. This rather important dependence 
on the parton densities can be considered as an extra handle provided by the asymmetry 
measurements in constraining the parton distribution functions. Conversely, with a more 
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LO 
NLO 


a(pb) 




A FB 


SA FB 




LHC 


o y i = 53 
57 

a Vb = 1.8 
1.9 


^FB 

A j 

r d 

^FB 


8.7 x 10~ 3 

6.8 x 10~ 3 
1.2 x 10~ 2 
1.1 x 10~ 2 
7.5 x 10~ 2 

4.9 x 10~ 2 


4.4 x 10~ 4 
4.2 x 10~ 4 
4.4 x 10~ 4 

4.2 x 10~ 4 

2.3 x 10~ 3 
2.3 x 10~ 3 


1.3 x 10~ 3 

1.3 x 10~ 3 
8.8 x 10~ 4 
1.1 x 10~ 3 
8.7 x 10~ 4 

1.4 x 10~ 3 


Tevatron 


a Vj = 8.0 
7.5 

a Vb = 0.04 
0.07 


4 cs 

^FB 

A j 

r d 

^FB 


6.4 x 10~ 2 

5.5 x 10~ 2 
9.9 x 10~ 3 
1.1 x 10~ 2 
5.5 x 10~ 2 
2.7 x 10~ 2 


3.5 x 10~ 3 

3.6 x 10~ 3 

3.5 x 10~ 3 

3.6 x 10~ 3 
5.1 x 10~ 2 

3.7 x 10~ 2 


1.4 x 10~ 3 
1.7 x 10~ 3 

8.1 x 10~ 3 

7.2 x 10~ 3 

2.5 x 10~ 2 
4.7 x 10~ 2 



TABLE III: Estimates for the e~e + j and e~e + b cross sections and asymmetries denned in the 
text with M e - e + in the range [75, 105] GeV. The first row of each entry is the LO result, while 
the second one refers to the NLO. The integrated luminosity as well as the cuts can be found in 
the text. The statistical precisions are also given, to be compared with the current effective weak 
mixing angle uncertainties at LEP and SLD from asymmetries only 1.6 x 10 -4 , and from Ap B at 
the Z pole 2.9 x 10" 4 [im- 
precise knowledge of them, the charge asymmetries can be used for precision measurements. 



CONCLUSIONS 



In summary, the large Vj production cross section at hadron colliders and the possibility 
of measuring the lepton asymmetries relative to the final jet allow for a precise determination 
of the effective electroweak mixing angle. We have evaluated them to NLO, confirming to a 
large extent the LO results. If there is an efficient 6-tagging and charge identification, these 
events with a b jet also allow for a new determination of Ap B . The corresponding statistical 
precisions are collected in Table IIHI As in Drell-Yan production [23| , these processes can be 
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also sensitive to new physics for large M e - e +, especially to new gauge bosons. 
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